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ABSmCT 


A cotuniter-oriented  method  for  obtaining  dynamical  equations  of 
motion  for  large  mechanical  systems  or  "chain  systems"  is  presented. 

A chain  system  is  defined  as  an  arbitrarily  assembled  set  of  rigid 
bodies  such  that  adjoining  bodies  have  at  least  one  common  point  and 
such  that  closed  loops  are  not  formed.  The  equations  of  motion  are 
developed  through  the  use  of  Lagrange's  form  of  d'Alembert's  principle. 

The  method  is  illustrated  and  applied  with  human-body  models 
and  finite-segment  cable  models.  The  human-bocfy  models  are  configured 
to  simulate  a crash-victim.  Results  with  several  applied  deceleration 
profiles  agree  very  well  with  available  experimental  data.  The 
cable  model  is  configured  to  simulate  an  off-shore  oil  rig  or  ship's 
crane  with  a partially  submerged  towing  cable. 
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I.  INIRDDUCTION 


Harry  mechanical  systans  and  devices  can  be  effectively  modelled 
by  systans  of  rigid  bodies.  If  a rigid-body  model  of  a mechanical 
systan  (called  a "finite-segment"  model)  consists  of  a system  of 
connected  (that  is,  non-dis joint)  linked  rigid  bodies  whidi  do  not  form 
closed  loops  or  circuits,  it  is  called  a "general  chain  system"  (or 
"open-chain")  Figure  1.  depicts  such  a svstem.  The  bodies  are  "linked" 
such  that  adjacent  bodies  share  at  least  one  coctmon  point,  thus  allow- 
ing for  either  hinge  or  ball-and-socket  connections.  Examples  of 
general  chain  svs terns  are:  hanan-bodv  models;  chain  and  finite- 
sepnent  cable  models;  manipulators;  and  finite-segiient  antenna  and 
beam  models. 

Recently  there  have  been  a number  of  attaipts  to  develop  efficient 
methods  for  obtaining  equations  of  motion  for  such  systems  [1-13]*. 

These  efforts  generally  proceed  by  first  modelling  or  replacing  the 
given  d5mamical  system  by  a discrete  system  or  chain  of  interconnected 
rigid  bodies.  n3mamical  equations  of  motion  are  then  written  for  the 
chain  system.  In  the  derivation  of  these  equations,  some  investigators 
use  Lagrange's  equations,  some  use  Newton's  laws,  and  some  use  other 
modified  geometrical  theoi  es.  But  each  has  the  objective  of 
efficient  derivation  of  computer-oriented  equations . The  relative 
advantages  and  disadvantages  of  these  various  methods  depends  unon: 


Ninfcers  in  brackets  refer  to  references  at  the  end  of  the  Report. 
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(i)  the  particular  dynamical  principle  which  is  used,  and  (ii)  the 
mediod  of  organizing  the  complex  geocnetry.  The  difficulties 
encountered  in  these  approaches  usually  include  some,  and  sometimes 
all,  of  the  follcwing:  (i)  the  introduction  of  non-vrorking  constraint 
forces  between  adjoining  bodies  (eq.  Newton's  laws);  (ii)  the  tedious, 
often  unwieldy,  calculation  of  derivatives  (eq.  Lagrange's  equations); 
(iii)  the  complex  geometrical  description  of  the  system;  and  (iv) 
the  solution  of  the  developed  eoviations. 

In  this  report  a method  of  obtaining  eouations  of  motion  is 
presented  which  systematically  avoids  each  of  these  difficulties . 

The  method  uses  Lagrange's  form  of  d'Alenbert's  principle  [1A-I7j 
which  provides  for  the  automatic  elimination  of  the  non-working 
internal  constraint  forces  without  introducing  tedious  differentiation 
or  other  calculation.  The  method  also  uses  a geometrical  organization 
and  accounting  procedure  as  developed  by  Kane  [16]  and  Huston  and 
Passerello  [8,10,12,18,19].  The  method  allcws  the  system  to  be  in 
any  general  force  field  and  either  the  moments  or  the  orientation 
between  adjoining  bodies  may  be  specified  or  left  unknown.  Finally, 
the  method  leads  to  governing  equations  which  may  easily  be  applied 
with  any  general  chain  system  sxach  as  human-body  models , finite- 
segment  cable  models,  manipulor  models  or  flexible  beam  models. 
Furthermore,  the  form  of  the  governing  equations  is  ideally  suited 
for  developing  computer  algorithms  for  obtaining  numerical  solutions. 

The  report  basically  outlines  the  results  of  research  under  the 
support  of  the  Office  of  Naval  Research  Contract  N00014-75-C-1164.  It 
contains  restilts  of  the  application  of  the  abcrve  method  with  human-body 
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models  in  high-acceleration  environments  (crashes)  and  with  partially 
submerged  tewing  cables. 


The  balance  of  the  report  is  devided  into  six  parts  with  the 
following  part  presenting  some  geometrical  background,  notation,  and 
other  preliminary  ideas  useful  in  the  analysis . The  kinematics , 
dynamics,  and  governing  equations  are  developed  in  the  next  two 
parts.  The  application  with  himan-bodv  models  and  cable  models  are 
presented  in  the  subsequent  two  parts,  and  a brief  discussion  and  set 
of  conclusions  are  presented  in  the  final  part. 


II.  PRELIMINARY  CONSEDERATIO^ 


Body  Connection  Array 

Consider  t±ie  chain  system  shown  in  Figure  1.  To  organize  the 
gecxnetrical  accounting  of  this  system,  select  one  body  of  the  system 
as  the  reference  body  and  call  it  Next,  nutrber  or  label  the  other 
bodies  of  the  system  in  ascending  progression  away  from  B^  as  shown 
in  Figure  2.  The  configuration  and  kinematics  of  each  body  of  the 
system  may  now  be  developed  relative  to  B^  which  in  turn  has  its 
configuration  defined  relative  to  an  inertial  reference  frame  R 
(Figure  2.)*. 

Althou^  this  nunhering  scheme  does  not  lead  to  a unique 
labelling  of  the  bodies,  it  can  nevertheless  be  used  to  describe  the 
chain  structure  throxjgh  the  "body  connection  array"  as  follows : 

Let  LCk) , k=l,  N be  an  array  listing  the  indices  of  the 

adjoining  lower-nunbered  body  for  each  body  Bj^.  For  example  for  the 
system  shown  in  Figure  2. , LCk)  is 

L(k)  = (0,1,2,2,A,1,6,1,8)  (2.1) 

where 

(k)  = (1,2, 3, 4,5, 6, 7,8, 9)  (2.2) 

and  where  0 refers  to  the  inertial  reference  frame  R.  It  is  not 
difficult  to  see  that  given  L(k)  one  could  readily  define  the  tocology 
or  arrangement  of  the  system.  That  is.  Figure  2.  could  be  constructed 
by  simply  knowing  L(k).  It  is  shown  in  Part  III  of  the  report  that  L(k) 


An  alternative  approach  is  to  reference  each  body  independently  to  R, 
but  this  is  found  to  be  very  inconvenient  in  describing  the  configuration 
of  actual  systems. 


can  be  useful  in  t±ie  development  of  expressions  of  kinematical 
quantities  needed  for  an  analysis  of  the  system's  dynamics. 


Shifter  Transformation  ffetrices 


Next,  consider  a typical  pair  of  adjoining  bodies  such  as  Bj  and 
as  shown  in  Figure  3.  The  orientation  of  Bj^  relative  to  Bj  mav  be 
defined  in  terms  of  the  relative  inclinations  of  the  dextral 
orthogonal  ijnit  vector  sets,  and  (i=l,2,3)  fixed  in  B^  and  Bj^ 
as  shown  in  Figure  3.  Specifically,  let  B.  and  B^^  be  oriented  such 
that  the  and  the  are  respectively  parallel.  Then  Bj^  may  be 
brought  into  any  orientation  relative  to  B^  by  three  successive  dextral 
rotations  about  axes  parallel  to  nj^2>  115,3  through  the  angles 

and  and  are  then  related  to  each  other  as 


q . . = SJK.  Hi 
ij  1 im  ikm 


(2.3) 


where  SJK  is  a 3 x 3 orthogonal  transformation  matrix  called  a 


'shifter”  and  defined  as  [20,21,22] 

SJK.  = n. . • q, 
im  1 l:<m 


(2,4) 


(Regarding  notation,  repeated  indices,  such  as  m in  the  ri^t  side 
of  Eouation  (2.3)  renresent  a sum  over  the  range  (eg.  1,  .,,,3)  of  that 


index.) 


SJK  may  itself  be  written  as  the  product  of  three  orthogonal 


transformation  matrices  as  [19] 
SJK  = oJK  SJK  yJK 
where  ctJK,  SJK,  and  y>JK  are 


(2.5) 
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1 


0 


0 


oJK  = 


0 -Soj^ 

0 Soj^  Cak 


0 

m=  0 10 

-s^  0 


(2.6) 


% -Sy^  0 
Sy^  Cy^  0 
0 0 1 

where  S and  C are  abbreviations  for  sine  and  cosine  respectively. 

From  Equation  (2.3)  it  can  be  seen  that  the  shifter  transformation 
matrix  obeys  the  following  identity  rules: 

SJJ  = I = SJK  SKJ  = SJK  SJK"^  (2.7) 

where  I is  the  identity  matrix.  Also , with  three  bodies , say  , 
and  the  shifter  obeys  the  chain  rule: 

SJL  = SJK  SKL  (2.8) 

These  expressions  may  be  used  to  transfer  conponents  of  vectors 
referred  unit  vectors  in  one  body  of  the  system  into  components  of 
vectors  referred  to  unit  vectors  in  any  other  body  of  the  system  and, 
in  particular,  to  the  inertial  reference  frame  R.  For  example,  if  a 

8 
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(2.9) 


typical  vectcfr  v is  expressed  as 

V = rii,-  “ ■ 

- 1 Iki  1 loi 


then 

= SOKy  (2.10) 

where  the  index  o refers  to  the  inertial  reference  frame  R.  (There 
is  no  sum  over  repeated  indices  in  parentheses.) 


Shifter  Derivatives 

Finally,  since  the  shifter  transformation  matrices  play  a 
central  role  througjiout  the  analysis,  it  is  helpful  to  also  have 
algorithms  for  their  derivatives,  especially  the  derivative  of  SOK. 
Such  an  algorithm  may  be  obtained  through  Equation  (2.4).  For  SOK, 
Ecruation  (2.4)  may  be  written  in  the  form: 


SOK. . = n . 
rj  loi 


Ikj 


(2.11) 


where  the  are  fixed,  and  therefore  constant,  in  R.  Hence, 
differentiating  in  Equation  (2.11)  leads  to 


d(SOK.  .)/dt  = n^. 


R 


dn^j  /dt 


(2.12) 


^diere  the  index  R in  a /dt  indicates  that  the  derivative  of 
n.  . is  conputed  in  R.  (See  Reference  [23].)  However,  since  the  n,  . 

~K.j 

are  fixed  in  their  derivatives  may  be  written  as  ^ where 
is  the  angular  velocity  of  Bj,.  in  R [23].  Equation  (2.12)  may  then 
be  written  as 


d(S0Ky)/dt  - 


(jj,  X n,  . = -(jj, 


X n . 

~oi 


n,  . 

~kj 


(2.13) 


or  as 


d(SOK)/dt  = WDR  SOK 
where  MDK  is  a matrix  defined  as 


SOK. 

im 


e. 

um 


(2.14) 

(2.15) 


and  where  are  the  conponents  of  referred  to  and  is  the 

standard  permutation  syirbol  [20,22] . (WDK  is  simply  the  matrix  whose 
dual  vector  [22]  is  !^.) 

Equation  (2.14)  shows  that  the  derivative  of  the  shifter  matrix 
may  be  computed  by  a matrix  multiplication  and  thus  is  ideally  suited 
for  development  into  a numerical,  computer  algorithm. 
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III.  KINEMAnCS 


Degrees  of  Freedom  and  Generalized  Coordinates 

The  chain  system  shown  in  Figure  1.  will,  in  general,  have 
3N  + 3 degrees  of  freedom.  These  may  be  defined  in  terms  of  3N  + 3 
generalized  coordinates  Xj  ....  3N+-3).  Let  X2»  arid  X3 

represent  the  position  coordinates  of  an  arbitrarily  selected 
reference  point  say  0^  of  in  R (See  Figure  4.) . Then  let  the 
remaining  3N  coordinates  be  divided  into  N triplets  of  coordinates 
representing  the  relative  orientation  angles  of  the  adjoining  bodies 
as  described  above.  For  example,  X41  X5.  and  Xg  define  the  orientation 
of  in  R and  X3;^+]_>  X3j^^2'  ^3k+3  define  the  orientation  of  B^^ 

relative  to  Bj  where  Bj  is  the  adjacent  lower  numbered  body  to  a 
t3T)ical  body  B^^  (See  Figure  3.) 

Angular  Velocities 

The  angular  velocity  of  a typical  body  of  the  system,  say  Bj^, 
relative  to  R is  readily  obtained  by  the  familiar  addition  formula 
[16,23]  as 

^,  = Si  + . . . ^ (3.1) 

where  the  terms  on  the  right  represent  angular  velocities  of  the 
subscripted  body  relative  to  its  adjacent,  lower-nurbered  body,  and 
xdiere  the  sun  is  taken  over  the  bodies  in  the  chain  from  Bi  outward 
through  the  branch  containing  Bj^.  For  exanple,  referring  to  Figure  2., 
lOc  is 

0)^  = 1^1  + ^2  + + 01^  (3.2) 
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\diere  is  the  angular  velocity  of  relative  to  R,  0^2  is  the 
angular  velocity  of  B2  relative  to  B^,  ^ is  the  angular  velocity  of 
B^  relative  to  B2,  and  is  the  angular  velocity  of  B^  relative  to  B^. 
(Regarding  notation,  the  "hats"  are  ijsed  to  designate  relative 
angular  velocities  with  respect  to  the  adjacent,  lower -nurbered  body, 
and  the  terms  without  the  hats  designate  absolute  angular  velocities 
(that  is,  relative  to  R) . Hence,  and  are  the  same.) 

Ihe  L(k)  array  introduced  above  can  be  ijsed  to  form  a 
convenient  expression  for  the  sim  in  Equation  (3 . 1) . To  see  this , 
consider  the  example  of  Equation  (3.2).  The  sijbscripts  on  the  right 
side  of  Equation  (3.2)  (that  is,  1,2, 4, 5)  may  be  obtained  from  L(k) 
as  follows:  Consider  L(k)  as  a function  mapping  the  (k)  array  into 
the  L(k)  array  (See  Eqijations  (2.1)  and  (2.2)).  Then,  using  the 
notation  that  L°(k)  = (k) , L^(k)  = L(k) , L^(k)  = La(k)),  L^(k)  = La^(k)), 
etc.,  it  is  seen  from  Equation  (2.1)  that 

L°(5)=5,  L^(5)=4,  L^(5)=2,  L^(5)=l  (3.3) 

Therefore,  may  be  written  as 

3 

*li5  ~ ^ where  q = L^(5)  (3.4) 

p=0 

Hence,  in  general,  the  angular  velocity  of  B^^  may  be  written  as 
r 

~ ^ where  q = (k)  (3.5) 

p=0 

and  where  r is  the  index  such  that  L^(k)=l. 
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Now,  in  view  Figure  3. , and  t±ie  description  of  the  relative 
orientation  of  with  respect  to  ^ nay  be  written  as  [19,21] 


?k  = ^Kl  + \^<2  + (3.6) 

where  and  ^ are  unit  vectors  parallel  to  '^k2> 

during  the  successive  dextral  rotations  defining  Oj^,  and 
as  described  in  the  foregoing  part  of  the  report.  In  tentis  of 
(i=l,2,3),  the  unit  vectors  fixed  in  B j , the  adjacent  lower 
nurrbered  body,  ^ (i=l,2,3)  are 

5cl  ~ !Ijl  ~ “W 


!?k2  ■ "~jm 

5k3  ■ 


(3.7) 


Hence,  ^ becomes 

or  in  terms  of  the  unit  vectors  fixed  in  R 

5k  = ( \^ml  + 


(3.8) 


(3.9) 


By  s\±)stituting  from  Equation  (3.9)  into  Eqijation  (3.5),  (Oj^  may 


be  written  in  the  form 


‘%i  ’^om 


(3.10) 


where  there  is  a sxm  from  1 to  3Nf3  on  2,  and  from  1 to  3 on  m.  From 


Equation  (3.9),  it  is  seen  that  the  non- zero  take  one  of  the 
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t±iree  forms 


= ^nn'^n2 


(3.11) 


depending  on  whetiier  is  t±ie  first,  second,  or  third  dextral  angle 

defining  the  orientation  of  relative  to  Bj . Also , from  Equation 

(3.5)  it  is  seen  (±iat  for  two  bodies  B_  and  B , in  the  same  branch  of 

r s 

the  chain,  for  r greater  than  s,  if  # 0. 


Lar  Accelerations 


The  angular  acceleration  of  Bj^  relative  to  R may  be  obtained  by 
differentiating  the  angular  velocity  expression  of  Equation  (3.10). 


Noting  that  the  n are  fixed  in  R,  this  leads  to 

^ -om  ’ 


(3.12) 


Frctn  Equation  (3.11),  the  non-zero  ^re  found  to  take  one  of  the 


three  forms: 


S^nn  “>^3 


(3.13) 


+ SOJ^  oJK  + 6JK„, 
tm  np  pj 


depending  on  whether  is  the  first,  second,  or  third  dextral  angle 
defining  t±ie  orientation  of  Bj^  relative  to  Bj . The  SQJ  are  given  by 
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•»<-' 


Equation  (2.14)  and  the  aJK  and  8JK  are  obtained  by  differentiating 
the  expressions  in  Equation  (2.6). 


Mass  Center  Velocities 

The  velocity  and  acceleration,  relative  to  R,  of  the  mass  centers 

of  the  bodies  of  the  system  may  be  obtained  as  follows:  First, 

recall  that  0^^  is  an  arbitrarily  selected  reference  point  of 

(See  Figure  4.,  lor  exanple) . Next,  let  be  the  connection  point 

or  cctnnon  point  of  two  typical  bodies,  say  B.  and  B.  where  B.  is  the 

J J 

adjacent  lower-numbered  body  of  Bj^  (k=2,  , . , , N) , and  let  Oj^  be 
called  the  "reference  point"  of  Bj^.  Then  let  be  the  position 
vector  of  relative  to  0^ . Finally,  let  be  the  mass  center  of 

Bj^  (k=l,  N),  and  let  r^  be  the  position  vector  of  relative  to 

Oj^.  (Ci^  is  thus  fixed  in  B^  and  r^^  is  fixed  in  Bj^.)  Hence,  the 
Txjsition  vector  of  Gj^  relative  to  a fixed  point  0,  in  R mav  be 
written  as 

u-1 

(3.14) 


Ek  - «oi »*Ih  I'kh  + I SOSth  W 'Joi 

vhere  s=L^(k),  S=L^^(k),  and  u is  the  index  such  that  L^(k)=l,  and 
v^iere  is  the  position  vector  of  0^^  relative  to  0.  Therefore,  bv 
differentiating  in  R,  the  velocity  of  Gj^  in  R _s 


u-1 


- (lot  + SdlKjj,  + Z *lhWEoi 

0=0 


(3.15) 


By  using  Equations  (2.14),  (2.15),  and  (3.10),  Vj^  may  be  written 
in  the  form 


YV.  ^kim  Yooi 


(3.16) 
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are 


IV.  EQUATIONS  OF  MDTION 


As  mentioned  earlier,  the  governing  d3matnical  equations  of 
motion  of  a general  diain  system  such  as  is  shown  in  Figure  1 . , 
can  be  obtained  conveniently  using  Lagrange's  Form  of  d'Alembert's 
principle  [14,17].  Since  the  non-working  internal  constraint 
forces  acting  between  the  bodies  of  the  system  are  automatically 
eliminated  in  using  this  principle,  only  the  externally  applied,  or 
active  forces,  acting  on  the  system,  and  the  so-called  inertia 
forces  of  the  system  need  be  considered  in  the  analysis. 

Kinetics 

Imagine  the  system  of  Figure  1.  to  be  subjected  to  an  arbitrary 
external  force  field.  Ihen  let  the  ensuing  forces  acting  on  a 
typical  body  of  the  system  be  replaced  bv  an  equivalent  force 
system  consisting  of  a single  force  passing  throu^  together 
with  a couple  with  torque  1^.  Then  the  so-called  "generalized 
active  force"  F^,  corresponding  to  the  generalized  coordinate  x^. 
acting  on  is  [16] 

■ ''kto.  '■km  + k m ^ ( -1 3Nf3;  k=l N)  (4.1) 

where  there  is  no  svm  on  k,  but  there  is  a sum  from  1 to  3 on  m. 

Fj^  and  are  the  components  of  ^ and  respectively.  Now, 

if  the  generalized  active  forces  of  each  of  the  bodies  of  the  system 
are  added  together,  the  resiilt  is  the  total  generalized  active  force 
for  the  entire  system.  Hence,  if  there  is  a sum  over  k from  1 co  N 
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in  Equation  (4.1),  then  represents  the  total  generalized  active 
force  on  the  system  for  the  generalized  coordinate  • • • ,3N+-3) . 

In  a similar  manner,  let  the  inertia  forces  on  a typical  body 
of  the  system  be  replaced  by  a single  force  passing  through 
together  with  a couple  with  torque  Ihen  ^ and  may  be  written 
as  [16]: 


and 


^ sun) 


Mk  = -^k  * ^ 


(4.2) 


(4.3) 


where  m^^  is  the  mass  of  and  is  the  inertia  dyadic  of  B^  relative 
to  (k=l,...,  N).  Ihrou^  use  of  the  shifter  transformation 


matrices  may  be  written  in  the  form 


5k  *”  ^nn  1cm  Ion 


(4.4) 


The  so-called  "generalized  inertia  force"  corresponding  to  the 
generalized  coordinate  x^,  acting  on  is  then  [16] 

V - "kM  ''kS  “Un  «kS  «-l k-l N)  (4.5) 

vdiere  there  Ls  no  sum  on  k,  there  is  a srm  from  1 to  3 on  m,  and  F^ 
and  are  the  components  o'*^  F^^*  and  respectively.  As  above, 
if  the  generalized  inertia  forces  of  each  of  the  bodies  of  the  system 
are  added  together,  the  resiolt  is  the  total  generalized  inertia 
force  for  the  entire  system.  Hence,  if  there  is  a sum  over  k from 
1 to  N in  Equation  (4.2),  then  F^^*  represents  the  total  generalized 
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inertia  force  on  the  system  for  the  generalized  coordinate 


t 


Governing  Equations 

Lagrange's  form  of  d'Alembert's  principle  states  that  the  sum 
of  the  total  generalized  active  force  and  the  total  generalized 
inertia  force,  for  each  generalized  coordinate  of  the  system,  is  zero. 
Hence,  the  governing  dynamical  equations  of  motion  for  the  system  are 
Fj^  + F^*  = 0 (il-l,...,3N4-3)  (4.6) 

By  STjbstituting  Equations  (3.10),  (3.12),  and  (3.20)  into 
Equations  (4.2)  and  (4.3)  and  ultimately  into  Eauations  (4.6),  the 
governing  equations  of  motion  may  be  written  in  the  form 

= a=l,....3Nf3)  (4.7) 

where  there  is  a sum  frcm  1 to  3Nf3  on  p and  where  a,^^  and  f^^^  are 

given  by 

~ '\^kpm\£m  \nn  ‘^^kpm  “^n  (^-8) 

and 

^kqm  ^q  ^kmn  ^Im  ^qn  ^q 

(4.9) 


®nnh  “^qn  ‘^sr  ‘^h  ^q  ^s^ 


where  there  is  a sm  from  1 to  N on  k,  from  1 to  3Nf3  on  q and  s,  and 
from  1 to  3 on  the  other  repeated  indices. 

Equations  (4.7)  form  a set  of  31^3  simultaneous  ordinary 
differential  equations  for  the  3Nf3  generalized  coordinates  Xj^  of  th® 
system.  If  some  (or  all)  of  the  Xj^  sere  specified,  the  differential 
equations  become  algebraic  equations  for  the  unknown  forces  or  moments 


20 


In  t±ie  following  two  parts  of  the  report  these  equations  are 
developed  and  solved  for  a htman-body  model  in  a variety  of  high- 
acceleration  configurations,  and  for  a finite-segment  cable  model. 
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V,  APPLICmON  WITH  HIJMAN-BODY  MDDELS 


As  mentioned  earlier,  a principal  application  of  the  foregoing 
analysis  is  studying  the  d3mainics  of  a human-body  model.  Indeed,  a 
desire  to  obtain  insight  into  the  d5maniics  of  space-walking  astronauts, 
athletes,  and  crash- victims  has  stimulated  much  of  the  development  of 
the  foregoing  analysis. 

Gross  motion  simulation  of  the  human  bodv  naturally  leads  to 
finite-segment  or  chain-systan  modelling  if  one  thinks  of  the  human 
body  in  terms  of  its  skeletal  structure.  That  is,  by  considering  the 
hands,  feet,  arms,  legs,  head,  and  torso  as  rigid  bodies  and  the 
muscles  and  ligaments  as  springs  and  dampers  acting  at  the  joints, 
the  ensxxLng  model  is  precisely  a general  chain  system  as  studied  in 
the  foregoing  parts  of  the  report.  In  this  part  of  the  report  a 
stmnary  of  research  resiolts  obtained  by  using  tlie  analvsis  in  the 
develooment  and  apnlication  of  a crash-victim  comouter  simulation 
code  [24-29],  comnonlv  known  as  "UCIN”,  is  presented. 

There  have,  of  course,  been  numerous  earlier  attempts  to  model 
the  d3mamics  of  the  human  body  and  in  particular,  the  d3mianiics  of  a 
crash  victim.  The  number  of  these  efforts  significantly  increased 
during  the  past  decade  with  the  availability  of  hi^-speed  digital 
computers.  Indeed,  there  even  exists  a number  of  recent  survey  papers 
[30-33]  outlining  this  work.  From  these  papers  it  appears  that  much 
of  the  significant  work  in  general  human  body  dynamics  and  modelling 
may  be  found  in  References  [18,19,34-38]  and  in  gross -motion, 
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crash-victim  simulation  in  References  124-33,39-78],  The  specific 
approaches  discussed  in  References  [18,19,24,29]  have  led  to  the 
development  of  the  foregoing  general  analysis  and  ultimately  to  the 
development  of  the  UCIN  model  as  described  in  the  following  section. 

The  UCIN  Mxiel 

The  model  consists  of  12  rigid  bodies  representing  the  human 
linbs  together  with  a vehicle  cockpit  as  shown  in  Figure  5.  The 
twelve  bodies  of  the  model  are  connected  together  with  ball-and-socket 
joints  vhich  allow  for  the  inclusion  of  springs  and  dampers  to 
simulate  the  human  connective  soft  tissue  svich  as  discs,  muscles,  and 
ligaments . 

Forty-five  variables  are  required  to  describe  the  position  and 
orientation  of  the  model.  These  are: 


X3_.X2.X3 


position  of  the  vehicle  relative  to 
an  inertial  frame 


X/^.X^.Xg 


orientation  of  the  vehicle  relative 
to  an  inertial  frame 


X7.X3.X9 


Xio»Xii.Xi2 


Xi3.Xi4,Xi5 


Xi6.X17.Xi8 


X19.X20.X21 


X22.X23.X24 


position  of  a reference  point  in  82, 
the  lower  torso  relative  to  the  origin  of 
the  vehicle  frame 

orientation  of  82,  the  lower  torso, 
relative  to  the  vehicle  frame 

orientation  of  8^,  the  middle  torso, 
relative  to  82,  me  Icwer  torso 

orientation  of  8, , the  upper  torso, 
relative  to  8^.  tne  middle  torso 

orientation  of  8c,  the  upper  left  arm, 
relative  to  Che  upper  torso 

orientation  of  8, , the  lower  left  arm 
relative  to  8^,  the  upper  left  arm 
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X25.X26.X27 


of  By,  t±ie  upper  right  arm, 

B^,  the  upper  torso 

of  Bo,  the  lower  right  arm. 

By,  the  upper  right  arm 

of  Bq,  the  head 
B^,  the  upper  torso 

of  ByQ,  the  upper  left  leg. 

By,  Cne  lower  torso 

of  B,^ , the  lower  left  leg, 

B^g.'^the  i^per  left  leg 

of  B^,  the  upper  right  leg. 

By,  ffie  lower  torso 

of  B,,,  the  lower  ri^t  leg, 

Biy.'^the  upper  right  leg 

All  of  these  variables,  except  for  the  position  variables 
X]^,  X2.  X3>  X7,  Xg,  and  Xg  are  orientation  angles  generated  by 
successive  rotation  of  adjacent  bodies  relative  to  each  other  as 
described  in  Part  II  of  the  Report.  The  first  six  variables  define 
the  motion  of  the  cockpit  or  vehicle  frame,  B^,  relative  to  the 
inertial  frame,  R.  These  variables  nust  be  specified  or  given. 

Also,  variables  X22.  X24.  X28.  X30.  X37.  X39,  X43.  and  X45  are  usually 
specified  as  zero  to  simulate  hinge  joints  at  the  elbows  and  knees . 

The  remaining  31  variables  may  be  either  specified  or  left  as 
unknowns.  If  the  variables  are  specified  (e.g.  X]^g=0) , the  required 
moment  needed  to  maintain  that  specification  (e.q.  M^g)  is  determined. 

The  model  accepts  arbitrary  specification  of  external  forces  and 
moments  on  each  of  its  bodies.  These  forces  and  moments  are  repre- 
sented on  each  body  by  an  equivalent  force  systan  consisting  of  a 
single  force  passing  through  the  mass  center  of  the  body,  together 
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X28.X29.X30 

X3i,X32.X33 

X34,X35,X36 

X37,X38.X39 

X4o.X4i,X42 

X43.X44,X45 


orientation 
relative  to 

orientation 
relative  to 

orientation 
relative  to 

orientation 
relative  to 

orientation 
relative  to 

orientation 
relative  to 

orientation 
relative  to 


wit±i  a couple. 


The  model's  initial  position  is  generally  in  an  erect  sitting 
position  as  shown  in  Figures  6.  and  7.  In  this  configuration,  all 
the  body  coordinate  axes  and  t±ie  vehicle  frame  are  aligned.  Thus  in 
this  configuration,  all  the  orientation  angles  are  zero. 

The  model  has  a seat  which  is  modelled  by  springs  as  shown  in 
Figure  9.  There  are  seven  "one-way"  springs  which  mav  exert  forces 
on  the  model  with  the  points  of  contact  being  the  mass  centers  of 
bodies  2,  3,  4,  9,  10,  and  12.  ("One-way"  springs  exert  forces  only 
while  in  compression.)  One-vrav  viscous  damper  stops  are  used  to 
limit  the  seat  deflection.  The  force  F generated  by  a damper  stop 
is  of  the  form 


F = ^ ^ 
^ 0 


if 

if 


X 1 Xo 


(5.1) 


where  k is  an  arbitrary  constant,  x is  the  spring  deformation 
(ccmpression) , and  Xq  is  sn  arbitrary  spring  compression  limit. 

The  model  has  a floor  or  foot  rest  which  is  modelled  as  a linear 


spring. 

The  model  provides  for  the  use  of  up  to  ten  restraining  belts 
modelled  as  springs  attached  at  arbitrary  points  between  the  cockpit 
and  the  bodies  of  the  model. 

The  ball-and-socket  connection  joints  of  the  model  are  provided 
with  angle  stops,  modelled  bv  one-way  dampers,  to  simulate  motion 
constraints  of  the  human  linfcs.  An  angle  stop  generates  a moment  ^ 
between  the  bodies  of  the  form 
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CQ-C^(a-a^)a  if  a>a^ 

M = Cq  if  -a2<ot<a^  (5.2) 

Cq+Cj^  (01-02)  if  ce<a2 

v^ere  a is  a rotation  angle,  and  ct2  arbitrarily  specified 
maxirnLim  and  tiiinimLin  values  of  a,  and  Cq  and  are  arbitrarily  specified 
constants. 

Finallv,  the  rodel  rrrovides  ■'^or  the  use  of  twelve  intrusion 
surfaces  or  planes  to  sirnulate  the  cockpit  or  vehicle  interior. 

These  intrusion  s\xrfaces  are  as  follows: 

1)  Left  windshield  5)  Lower  left  door  9)  Firewall 

2)  Front  windshield  6)  Upper  right  door  10)  Top  dash 

3)  Right  windshield  7)  Lower  right  door  11)  Front  dash 

4)  Upper  left  door  8)  Roof  12)  Bottom  dash 

The  location  and  inclination  of  these  intrusion  surfaces  are  determined 
by  the  specification  of  the  position  of  a point  in  the  surface  together 
with  the  conponents  of  a vector  normal  to  the  siirface.  -v. 

The  UCIN-CRASH  Computer  Code:  Input/Output 

The  following  brief  paragraphs  provide  a general  description  of 
the  input  data  recruired  and  the  output  data  provided  by  the  computer 
code.  Additional  details  may  be  obtained  from  References  [27]  and  [28]. 
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The  input  consists  of  the  following: 


Physical  parameters:  These  are  the  masses,  inertia  dvadics,  mass 
center  positions,  connection  point  oositions,  and  orientation  angle 
limits,  for  each  of  the  12  bodies  of  the  model. 

Cockpit  geometry:  This  consists  of  a normal  vector  and  a location 
point  for  each  of  the  12  intrusion  surface  nlanes  listed  above. 

Also,  the  floor  position  and  a spring  constant  of  the  floor  model  are 
part  of  the  cockpit  specifications. 

Cockpit  motion:  The  cockpit  displacement  and  rotation  relative  to 
an  inertia  frame  required  as  input.  Typically,  it 

is  tiseful  to  express  this  in  terms  of  the  linear  and  angular 
acceleration  of  the  cockpit.  The  computer  program  is  written  so  that 
the  cockpit  acceleration  components  may  be  "read  in"  by  simply 
specifying  the  acceleration  profile.  (This  is,  in  effect,  a piecewise- 
linear  approximation  to  an  acceleration  curve.)  Six  (three  translation 
and  three  rotation)  acceleration  profiles  or  curves  mav  be  anp loved. 

Spring  and  damping  constants:  These  include  seat  constants,  restrain- 
ing belt  constants,  orientation  angle  constants,  and  neck  parameters. 
Also,  the  attach  points  of  the  restraining  belts  are  included  as  part 
of  this  data. 
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InjH conditions : This  includes  the  initial  values  of  the  unknown 
variables  and  their  derivatives.  Also,  the  external  forces  and 
moments  (if  anv)  which  are  applied  to  the  bodies  of  the  models  must 
be  specified. 

Integration  parameters:  This  consists  of  constants  recnrLred.  bv  the 
fourth-order,  Runge-Kutta  integration  technique  (RKGS)  and  it 
includes  the  starting  time,  the  ending  time,  the  step  size,  and  the 
error  permitted. 

The  output  consists  of  two  parts:  The  first  is  simply  a copy 
or  "echo"  of  the  input  data.  The  second  contains  at  each  output  step 
the  following: 

1)  The  value  of  all  variables  and  their  first  and  second 
derivatives . 

2)  The  joint  and  mass  center  positions  in  both  inertia 
space  and  relative  to  the  cockpit. 

3)  The  mass  center  velocities  and  accelerations. 

4)  The  mcments  and  forces  associated  with  variables 
which  are  specified, 

5)  Restraining  belt  forces, 

6)  Collisions  or  "hits"  with  intrusion  surfaces . 
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Validation  of  the  Model  and  Examples 


There  is  little  experimental  data  available  to  date  which  can 
be  used  to  check  or  verify  the  above  computer  code  and  others  like 
it.  However,  King,  Padgaonkar,  and  Mital  [79]  have  recently  con- 
ducted a series  of  experiments  with  a cadaver  in  an  impact  seat 
for  the  purpose  of  validating  the  conputer  model.  The  results  show 
remarkable  agreement  between  the  model  and  the  experiment  as 
shown  in  Reference  [79]. 

In  earlier  tests,  Begeman,  King,  and  Prasad  [80]  placed  a 
cadaver  in  an  impact  sled  and  they  measured  shoulder  belt,  vertical 
lap  belt  and  seat  pan  forces.  This  same  test  was  simulated  with 
the  UCIN  ccnputer  model.  For  the  deceleration  profile  shown  in 
Figure  8.,  the  experimental  and  computed  shoulder  belt,  vertical 
lap  belt , and  seat  nan  forces  are  shown  in  Figures  9 . , 10 . , and  11 . 

In  another  configuration,  experimental  data  from  a vehicle 
striking  a guard  rail  or  roadside  barrier  [81]  was  used  as  input  for 
the  computer  code.  For  the  specLfic  deceleration  profile  shown  in 
Figure  12.,  the  resu].tant  relative  displacanent  of  the  head  and 
chest  using  both  lap  belts  and  a combination  of  both  lap  and  shotilder 
felts,  is  shown  in  Figure  13. 

Finally,  in  an  attempt  to  measure  the  relative  effectiveness 
of  lap  and  shoxilder  belts,  a run  was  made  simulating  a front  end 
collision  of  a vehicle.  The  head  pitch  (forward  rotation)  of  the 
model  was  calculated  using  a lap  belt  and  a combination  lap  and 
shoulder  belt.  The  results  shown  in  Figure  14.  illustrate  the 
"whiplash"  effect  when  only  lap  belts  are  lised. 
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SLED  DECELERATION 


Figure  9.  Co»Tparison  of  Shoulder  Belt  Forces 


SHOULDER  BELT 


Figure  11.  Comparison  of  Se 


SEAT  FORCE 


Figure  12.  Vehicle  Deceleration 
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Figitre  13.  Copparison  Relative  Head  Potation  for  I.an  Pelts  and  Ror  a Combination  of  Lao  and  Shoulder  Belts. 
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Figure  14.  Corparison  of  I^ap  and  Shoulder  Belts  (Front  End  Collision) 
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HEAD  PITCH  RELATIVE  TO  UPPER  BODY 
LAP  AND  SHOULDER  BELTS 


VI.  APPLICATICN  wrm  FINTIE- SEGMENT  CABLE  MDDELS 


A second  application  area  of  the  foregoing  general  analysis 
is  the  study  of  the  dynamics  of  long,  heavy  cables.  As  with  human 
body  models,  a finite-segnent  cable  model  is  readily  identified  as 
a general  chain  system.  Indeed,  modelling  a cable  by  a finite- 
segment  model  sinply  involves  the  suibstituting  of  a linked  chain 
for  the  cable  as  shown  in  Figure  15.  This  model  has  the  obvious 
advantage  of  being  "linear",  that  is,  not  possessing  any  "branches". 
This  in  turn,  provides  a simolification  in  the  governing  eouations. 

Cable  dynamics  has  been  of  interest  to  researchers  for  some 
time.  But  recently,  with  the  advent  of  high-sneed  digital  computers 
and  finite-element  methods  and  finite-segment  modelling,  there  has 
been  increasing  interest  and  research  effort  in  cable  phenomena. 

Five  years  ago  Choo  and  Casarella  [82]  published  an  excellent  survey 
of  the  literature  and  analytical  methods  for  cable  dynamics.  They 
indicate  that  of  four  methods  for  studying  cable  dynamics  (the 
method  of  characteristics,  the  finite-element  method,  the 
linearization  method,  and  the  equivalent  lumped  mass  method) , the 
finite-element  method  is  the  most  versatile.  Indeed,  they  indicate 
that  the  finite-element  or  finite-segment  method  offers  the  best 
hope  for  a sinple  method  that  can  solve  nonlinear,  unsteady  state 
problems  with  good  accuracy  and  yet  require  only  a moderate  amount 
of  computation  time. 
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References  [83-92]  provide  a sutimary  of  recent  finite-element 
and  finite-segnent  approaches  to  cable  dynamics.  Ihe  approach  in 
this  report  is  to  make  a finite-segment  model  of  the  cable  as  in 
Figure  15.  and  to  then  follow  the  analysis  as  reported  in  Reference 
[93],.  Ihe  basic  d37namics  of  this  model  is  then  a specialization  of 
the  foregoing  general  theory. 


Ecpjations  of  Motion,  Computer  Code,  and  Numerical  Solutions 

The  governing  dynamical  eqi^ations  of  motion  are  of  the  exact 
same  form  as  those  presented  in  Part  IV  of  the  Peixjrt.  Indeed,  the 
only  modification  required  in  the  foregoing  analysis  is  a simplifi- 
cation in  the  form  of  the  mass  center  oosition  and  velocity  expressions 
(due  to  the  "linearity”  of  the  cable  model) , and  a specialization  of 
the  generalized  active  forces  to  account  for  the  fluid  drag  forces. 
Specifically,  the  expression  for  the  position  of  the  mass 


center  of  link  in  Equation  (3.14)  is  replaced  with 


(6.1) 


Then  the  velocity  of  the  mass  center  in  R (Eqijation  3.15)  becomes 


Ol 


(6.2) 


This  leads  to  the  following  expressions  for  the  non-zero 
(replacing  Equations  (3117),  (3.18)  and  (3.21)) 


kJJm 


'^£m  Ck=l,...,N;  )!.,iit=1,2,3) 


(6.3) 
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(6.4) 


k-1 


and 


(k=l,...,N;  2,=l,...,3N+-3;  tn=l,2,3) 


k-1 


^kJlm 


(6.5) 


(k=l,...,N;  ^=l,...,3Nf3;  tif=1,2,3) 

where  N is  t±ie  nurrber  of  links  or  segnents  in  the  model. 

Regarding  the  generalized  active  forces  due  to  the  fluid  drag, 
the  fluid  forces  on  the  cable  links  are  modelled  as  follows : If 
is  a typical  link,  the  fluid  forces  on  are  taken  to  be  equivalent 

to  a single  force  passing  through  (the  mass  center  of 

and  perpendiaiLar  to  the  axis  of  B^.  (Uie  axis  is  the  line  joining 
the  connection  points  of  \ and  Bj^,  .)  A sketch  of 

is  given  in  Figure  16.  Analytically  Fj^  is  expressed  as 


Ikw  = -p^^IyIclIYIcx 


(6.6) 


where  p is  the  fluid  mass  density,  Aj^  is  the  projected  link  area  on 
a plane  containing  the  axis  of  Bj^,  Cp  is  the  drag  coefficient,  and 

Yki  = ik  ='  ’^k' 


where  \y.  is  a unit  vector  parallel  to  the  axis  of 
velocity  of  Gj^,  and  v^  is  the  velocity  of  the  fl\rLd  (See  Figure  16.) 


vj^  is  the 


To  obtain  the  generalized  active  forces  associated  with  the 
fliiid  drag  (as  well  as  the  gravitational  forces)  is  substituted 
into  Equation  (4.1).  This,  in  turn,  leads  to  governing  dynamical 
equation  of  motion  of  the  form  of  Equations  (4,7) . 

As  with  human-body  crash-victim  ncdels , a cctnputer  code  has 
been  written  to  evaluate  the  coefficients  of  the  governing  equations . 
As  input,  the  code  requires:  the  nutiher  of  links;  the  masses; 
the  centroidal  principal  inertia  matrices;  the  mass  center  positions; 
the  connection  point  positions;  the  motion  profile  for  those  links 
with  specified  motion;  the  initial  configuration;  the  airbient  fluid 
velocity;  the  fluid  surface  hei^t  relative  to  R;  the  fltiid  mass 
density;  the  fluid  kineinatic  viscosity;  the  projected  link  areas; 
the  link  diameters ; and  the  mass  densities  of  the  links . The  code 
provides  for  the  evaluation  of  C^,  the  drag  coefficient  of 
Equation  (6.6)  from  an  algorithm  which  models  Hoemer's  [94]  drag 
coefficient  curve. 

The  governing  equations  of  motion  are  numerically  integrated 
using  a fourth  order  Runge-Kutta  technique.  The  output  of  the 
conputer  code  then  includes:  the  values  of  all  variables  and  their 
first  derivatives;  the  connection  point  positions;  the  mass  center 
positions ; the  mass  center  velocities  and  accelerations ; and 
the  momaots  and  forces  associated  with  the  specified  variables.  All 
of  the  output  is  given  at  arbitrarily  spaced  time  intervals. 
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Example  Application 

Consider  a rotating  surface  crane  with  a 25  ft.  boom  dragging  a 
50  ft.  cable  attached  to  a suhnerged  1 ft.  diameter  sphere,  as 
depicted  in  Figure  17.*  The  cable  is  modelled  bv  10  cvlindrical 
links,  each  5 ft.  long.  Ihe  cable  dianeter  is  1 in.  and  each 
link  has  a mass  of  0.4025  slug.  The  sphere  mass  is  7.727  slug. 

The  water  is  calm  and  its  surface  is  10  ft.  below  the  boom. 

The  boon  makes  a 90°  turn  in  5 sec.  The  angular  acceleration 
is  given  by  the  graph  in  Figure  18.  The  resultant  motion  of  the 
sphere  is  shown  in  Figures  19,  and  20.  To  determine  the  effect  of 
the  water  drag  force,  the  same  run  without  the  water.  It  is  seen 
(as  one  would  expect)  that  the  drag  forces  tend  to  counteract 
the  inertia  forces. 


Although  this  simulates  an  off-shore  oil  rig  or  a ship's  crane,  it 
does  not  represent  any  specific  physical  situation.  Indeed,  the 
intention  is  simply  to  demonstrate  the  kinds  of  problems  which  can 
be  studied  with  the  above  described  analysis  and  conputer  code. 
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Figure  18.  Angular  Acceleration  of  the  Boom 
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Flpjiire  19.  R/iclial  Position  of  the  Sphere 
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VII.  DISCUSSION  AND  CONCLUSIONS 


A new,  efficient,  computer-oriented  method  for  obtaining  ecuations 
of  motion  for  large  mechanical  systems  (chain  systems)  has  been 
developed.  These  equations  may  be  succinctly  written  in  the  form  of 
Equations  (4.7)  with  the  coefficients  determined  by  the  fovnr 
kinematical  matrices  and  These  matrices, 

in  turn  are  determined  by  simple  multiplication  algorithms  which 
enable  the  entire  procedure  to  be  reduced  to  a systematic  routine. 

The  development  of  this  method  was  made  possible  through  the 
use  of  Lagrange's  form  of  d'Alenfcert's  principle  which  provides  for 
the  automatic  elimination  of  nor-wrking  constraint  forces . 

Specifically,  it  is  this  feature,  combined  with  the  use  of  vector- 
matrix  notational  schemes,  which  make  possible  the  explicit  writing 
of  the  equations  of  motion.  Also,  the  analysis  makes  use  of  "local" 
as  opposed  to  "global"  coordinates.  That  is,  the  generalized 
coordinates  of  the  systan  are  (except  for  the  translation  of  the  first 
body  (or  segment)  relative  orientation  angles  measiired  between  the 
respective  bodies  as  opposed  to  absolute  orientation  angles  of  the 
bodies  (or  segments)  in  space.  This  allows  for  a more  convenient 
specification  of  initial  conditions  and  constraints,  and  for  a more 
convenient  interpretation  of  the  res\alts. 

The  range  of  application  of  the  analysis  and  procedures  is 
very  broad.  The  two  example  areas  of  h\jman-body  models  and  cable 
models  are  simply  two  areas  where  there  is  current  interest. 

Information  regarding  other  possible  areas  of  application,  the  use 
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of  program  tapes,  and  user  manuals  may  be  obtained  from  the  authors. 

Future  work  will  involve  extending  the  procedures  to  allcw  for 
translation  between  the  bodies  and  for  closed  loops  in  the  system 
models . 
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